Group properties and invariant solutions of a sixth-order thin film equation 

in viscous fluid 
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Using group theoretical methods, we analyze the generalization of a one-dimensional sixth-order 
thin film equation which arises in considering the motion of a thin film of viscous fiuid driven by 
an overlying elastic plate. The most general Lie group classification of point symmetries, its Lie 
algebra, and the equivalence group are obtained. Similar reductions are performed and invariant 
solutions are constructed. It is found that some similarity solutions are of great physical interest 
such as sink and source solutions, travelling-wave solutions, waiting-time solutions, and blow-up 
solutions. 
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I. INTRODUCTION 

In the past several decades there is an increasing inter- 
est in physics and mathematics literatures in higher-order 
nonlinear diffusion equations because they are models of 
various interesting phenomenon in fluid physics and have 
surprising mathematic structure and properties. Proba- 
bly, one of the most famous example is the fourth-order 
thin film equation in the form 



a > 



(1) 



which was first introduced by Greenspan in 1978 This 
equation describes the surface-tension-dominated motion 
of thin viscous films for the film height u{t, x) and spread- 
ing droplets in the lubrication approximation [J] . In par- 
ticular, for a = 3 it describes a classical thin film of 
Newtonian fluid, as reviewed in [2|, a = 1 occurs in the 
dynamics of a Hele-Shaw cell Q and a = 2 arises in a 
study of wetting films with a free contact line between 
film and substrate [1] . There also exist many interesting 
generalizations of the famous equation ([T]) (see @ and 
reference therein). 

Apart from the fourth-order equations, another inter- 
esting higher-order diffusion model is the sixth-order non- 
linear thin film equation in the form 



ut = {u" 



(2) 



which appear in flow modeling. The case m = 3, for in- 
stance, was first introduced by King in \Jj as a model of 
the oxidation of silicon in semiconductor devices or for 
a moving boundary giveri by a beam of negligible mass on 
a surface of a thin film [9] . Here u > will be treated as 
the thickness of a fluid film beneath an elastic plate and 



P — Uxxxx as the pressure within the film ^. The other 
derivatives of u can in the usual way be assigned different 
physical meaning, for instance T = —u^x is the bending 
moment on the overlying plate and S = Uxxx is the shear- 
ing force [l(j |. here all such expressions are dimensionless. 
An equation of this type can be used to model the mo- 
tion of a thin film of viscous fluid overlain by an elastic 
plate [8^; see also Hobart et al. [ll| and Huang et al. [l3] 
for possible applications of such modelling approaches to 
the wrinkling upon annealing of SiGe films bonded to Si 
substrates. Other plausible applications of Eq. ([2]), and 
suitable generalizations thereof, include a simple model 
for the influence of a crust on a solidifying melt or for a 
microfluidic pump (see Koch et al. [l^, for instance). 

Eqs. Ill and ([2]) are also the second and the third mem- 
ber of a hierarchy arising from the generalized Reynolds 
equation 



Ut = {u^Px)x, 



(3) 
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under different driving forces respectively. For gravity 
driven flows, we have p — u, giving the very widely stud- 
ied porous-medium equation (see, for example, Aronson 
[l^). For surface-tension driven flows we have p — —u^x, 
leading to the fourth-order thin film equation ([1]). For 
elastic plate driven flows, we have p = u^xxx, which give 
the sixth-order thin film equation ^ 8] . 

Up to now, the mathematic structure and proper- 
ties of the fourth-order thin film equation ([1]) have been 
widely investigated, including (non-) uniqueness, wetting 
behaviour and contact line motion, in particular opti- 
mal propagation rates, waiting time or dead core phe- 
nomena and self-similar solutions(see Hulshof [l^, for 
instance). Recent years, there are also many researches 
devoted to symmetry group structure and exact solu- 
tions of the fourth-order thin film equations ([T|) and their 
generalizations p^l - l25| . or searc hing for special invariant 
finite vector spaces of solutions [2g| . 

However, the sixth-order thin film equation ^ has 
been much less extensively investigated. It was only a few 
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researches that were devoted to quahtative mathematic 
properties such as the existence of weak solutions, ini- 
tial boundary value problems (see Bernis Friedman [27l |. 
King et al lisi, Smith et al. ^2^, Evans et al. Bar- 
rett et al. [31] for existing studies, the first two being 
analytical and the others primarily numerical), while the 
symmetry group properties and corresponding algebraic 
structure as well as explicit exact solutions of Eq. ([2]) still 
remain open. Therefore, the aim of the present work is 
to find such group properties, algebraic structure and ex- 
act solutions. To do this, we investigate alternatively a 
more general sixth-order nonlinear diffusion equation in 
the form 

Ut = {f{u)Uxxxxx)x, (4) 

than the original equation where f{u) is an arbi- 
trary smooth function depending on the geometry of the 
problem and /« ^ 0( i.e., is a nonlinear equation). 

We use the method of Lie groups, one of the powerful 
tools available to solve nonlinear PDEs, and which was 
discovered and applied firstly by S. Lie in the nineteenth 
century, but only in the last decades has it become a com- 
mon tool for both mathematicians and physicists (see for 
examples [s^-Hlj). The method consists of looking for the 
infinitesimal generators of a group of point transforma- 
tions which leave the equation under study invariant. An 
important point of the Lie theory is that the conditions 
for an equation to admit a group of transformations are 
represented by a set of linear equations, the so-called "de- 
termining equations" , which are usually completely solv- 
able. Having once found the groups of transformations, 
one can obtain a number of interesting results, which in- 
clude the possibility to reduce a partial differential equa- 
tion with two independent variables to an ordinary differ- 
ential equation with one independent variable, etc.. Solv- 
ing these reduced equations, one can obtain some partic- 
ular solutions for the original equations. These particular 
solutions are usually called "similarity solutions" or "in- 
variant solutions" p32H35| . When the equation contains 
"arbitrary elements" (a variable coefficient deriving from 
the particular equation of state chosen to characterize 
the physical mechanism. "Arbitrary elements" are func- 
tions or variable parameters whose form is not strictly 
fixed and can be assigned freely on the grounds of phys- 
ical hypotheses about the nature of the medium under 
consideration.), the theory gives rise to the problem of 
group classification of differential equations which is the 
core stone of modern group analysis (ssl . [35| . In par- 
ticular, in the past several years, a numbers of novel 
techniques, such as algebraic methods based on subgroup 
analysis of the equivalence group [42| - |45| . compatibility 
and direct integration [H, \^ (also referred as the Lie- 
Ovsiannikov method) as well as their generalizations (eg. 
method of furcate split [i^l, additional and conditional 
equivalence transformations fisi, |49| , extended and gen- 
eralized equivalence transformation group, gauging of ar- 
bitrary elements by equivalence transformations [lo, [HI) 
have been proposed to solve group classification problem 



for numerous nonlinear partial differential equations. Al- 
though a great deal of classification was solved by these 
methods, almost all of them are limited to the equations 
whose order are lower than four (see [5l| for details). 

In this paper we extend these new techniques, specific 
compatibility and direct integration as well as equiva- 
lence transformation techniques, to sixth-order nonlinear 
diffusion equations. We first carry out group classifica- 
tion of Eq. (j4l) under the usual equivalence group. The 
Lie group of point symmetries of Eq. , as a special case 
of Eq. and its Lie algebra are also obtained. Then 
similar reductions of the classification models are per- 
formed and invariant solutions are also constructed. It 
is found that some similarity solutions are solutions with 
physical interest: sink and source solutions, travelling- 
wave solutions, waiting-time solutions and blow-up solu- 
tions. 

The rest of this paper is organized as follows: In Sec.HIl 
we derive the equivalence group and perform the group 
classification related to Eq. (|4]). In Sec. IIIIl similar re- 
ductions of classification models are carried out. Sec. 
IIVI contains examples of some specific exact solutions, 
including sink and source solutions, travelling- wave solu- 
tions, waiting-time solutions and blow-up solutions, while 
in Sec. |V] some concluding remarks are reported. 



II. SYMMETRY CLASSIFICATION 

Background and procedures of the modern L ie group 
theory are well described in literature (32| - [35l l48l l49j |. 
Without going into the details of the theory, we present 
only the results below. 

Let 

Q = r(t, X, u)dt + £,{t, X, u)dx + (l)(t, x, u)du 

be a vector field or infinitesimal operator on the space 
of independent and dependent variables x, u. A local 
group of transformations G is a symmetry group of Eq. 
(g]) if and only if 

pr(6)Q|(A)=0, (5) 

whenever — ut~ \_f {u)uxxxxx\ = for every generator 
of G, where pr'^^^Q is the sixth-order prolongation of Q. 
Expanding Eq. ^ we get 

0* =(j)f"{u)UxUxxxxx + f'{u)UxxxxxV + f {u)Ux(jf '''''''' 
+ <iynu)Uxxxxxx + /(U)^""""" (6) 

which must be satisfied whenever Eq. @ is satisfied. 
Substituting the formulae of (j)^ , ^j^xxxxx ^xxxxxx 
into Eq. ([6]) we get an equation of t, x, u and the deriva- 
tives of r, ^, 0, u. Replacing ut by the right hand side of 
Eq. (HI) whenever it occurs, and equating the coefficients 
of the various independent monomials to zero, we obtain 
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the determining equations 

3(20,„-5e.,)/(u) + 0,/'H = 

{^xu-2^xx)f'{u)=Q 
'^4^xxu ^^xxx — 
{(l^xxu - ixxx)f'{u) = 

(rt-6ex)/H + 0/'H=O 

^4^xxxu ^Cxxxx — 
4>xxxxxxf{u) - 0t = 

4^xxxxx f ('^) ~t~ ~i~ ^^xxxxxuf ^xxxxxxf {'^^ — 

^4^xxxxu '^Cxxxxx — 
(^^^xxxxu ^xxxxx^ f (^) — 
^ (^4^xxxu ^xxxx^f ('^) 

The first three equations imply that (jj^u = £.xx — 0, 
which together with the sixth, the eighth and the ninth 
equations imply that 4>t = fpx — £,t = 0, so the determin- 
ing equations reduce to 



Tx = Tu = 
it — Cu — ixx — 
t>t^(l)x^ 4>uu = 

.(n-6e.)/(u) + 0/'(u) = O, 



(7) 



which is equivalent to 



< 



£^ — ax + h 
T — ct + d 
(j) = pu + q 
[ic-6a)fiu) + {pu + q)f'{u) = 



(8) 



where a, b, c, d, p, and q are arbitrary constants. 

In order to make the classification as simple as possible, 
we next look for equivalence transformations of class (j4]) , 
and then solve system ([S]) under these transformations. 
An equivalence transformation is a nondegenerate change 
of the variables i, x and u taking any equation of the 
form Q into an equation of the same form, generally 
speaking, with different f{u). The set of all equivalence 
transformations forms the equivalence group G"^ . To find 
the connected component of the unity of G"^ , we have to 
investigate Lie symmetries of the system that consists of 
Eq. Q and some additional conditions, i.e. 



fu'^x^xxxxx ^ f'^XXXXXXj 

/t = 0, 
=0. 



(9) 



That is to say we must seek for an operator of the Lie 
algebra A"" of G~ in the form 



X —T{t, X, u)dt + i{t, X, u)dx 

+ X, u)du + il>{t, X, u, f)df. 



(10) 



Here u and / are considered as different variables: u is on 
the space and / is on the extended space {t,x,u). 



The coordinates r, ^, of the operator (|T0)) are sought 
as functions of i, x, u while the coordinates ^ are sought 
as functions of i, x, u and /. 

Applying pr^^^X to Eq. © we get the infinitesimal 
criterion 



4^ — '^x'^xxxxx'4^ ~t- fu^xxxxx4^ 

+fuUxr^^^^ + u^xxxxx^ + fr^---' 



(11) 



which must be satisfied whenever Eq. (|9]) is satisfied 
Substituting the formulae of 0*, 0^, 0^^^^^^ ^xxxxxx^ 
-0^, and "0" into Eq. (fTT|) we get equations of t, x, u 
/, and the partial derivatives of t, ^, 0, u, /, and 
Replacing Ut, ft and /j. by the right hand side of Eq, 
([9|) whenever they occur, and equating the coefficients of 
various independent monomials to zero, we obtain 

f ixx =0, 6=0, = 

4>x = 0, (f>t ^ 0, 4>uu = 
= 0, iJt^ 0, = 

which can be reduced to 



(12) 



T = C^t + Cl 

C = C5X + C2 

(j) = ceu + C3 
^ip = (6C5 - C4)/ 

where ci, C2, . . . , ce are arbitrary constants. 
Thus the Lie algebra of G~ for class ([4]) is 

A~ = {dt,dx, du, tdt - fdf,xdx + 6/9/, 



Continuous equivalence transformations of class ^ are 
generated by the operators from In fact, G~ contains 
the following continuous transformations: 



t — t£4 + 61 , 

u = uee + £3, 



+£2, 



'5j 



where £1, £2, . . . , £6 are arbitrary constants. 

Solve the last equation of system dH) under the above 
equivalence group we can obtain three inequivalent 
equations of class (|4]) with respect to the transformations 
from G~: 

Case 1: f{u) is an arbitrary nonconstant function, 
the symmetry algebra of class Q is a three-dimensional 
Lie algebra which is generated by the operators 



Qi = dt, Q2 = dx, Q3 = tdt 



xdx. 



(13) 



Case 2: / = e^" mod G~{\ ^ 0), the symmetry al- 
gebra of class (jH) is a four-dimensional Lie algebra which 
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is generated by the operators 
Qi = jdu - tdt, 



B. 



f{u) = e^" (A / 0) 



(14) 



Case 3: / = u™ mod C{m ^ 0), the symmetry al- 
gebra of class (|4]) is a four-dimensional Lie algebra which 
is generated by the operators 



Qi = —udu - tdt, Q2 = dt, 

m 

Qz = -xdx udu, Qi = d^. 

m 



(15) 



From the above results, it is easy to see that equation 
^ is exactly corresponding to case 3, thus possess a 
four-dimensional symmetry algebra. 



III. SIMILARITY REDUCTION 

In order to obtain all the inequivalent reductions, we 
look for the one-dimensional optimal systems (see |33|). 
These systems, similarity variables and reduced equa- 
tions are listed below. In the following tables HI |IT1 IIIIl 
each row shows the infinitesimal generators Qi of each 
optimal system, as well as its similarity variable, simi- 
larity solution and reduced equation, a is an arbitrary 
constant, while A is a non- vanishing arbitrary constant. 
Note that in the case /(m) = u™ which corresponds to 
Eq. ([5]), we only consider m ^ 0, otherwise the equation 
is linear. 



A. /(it) is an arbitrary nonconstant function 

In this case, the symmetry operators are Eq. ([T^ . 
These operators satisfy the commutation relations 



[Qi, Q3] [Q2, Qz 



Q2 



and thus the corresponding symmetry algebra is a real- 
ization of the algebra ^3^5(0 < \a\ < 1) [52]. According 
to the results of Patera and Winternitz ^52] , an optimal 
system of one-dimensional subalgebras is those spanned 

by 

Qi, Q2, Qs, Qi + aQ2. 

Therefore, the corresponding similarity variables and re- 
duced ODEs can be easily calculated. Such results are 
Usted in Table HI 

Table I. Reduced ODEs for arbitrary nonconstant f{u) 

(let E = [f(v}Vyyyyy]y). 



i 


Subalgebra 


Ansatz u = 


y 


Reduced ODE; 


1 


(Qi> 


■"(y) 


X 


E = 


2 


(Q2> 


•"{y) 


t 


Vy =0 


3 


(Q3> 


■"{y) 




E = -yvy/6 


4 


(Qi +aQ2) 


■"(y) 


X — at 


E = —avy 



In this case, the symmetry operators are given by Eq. 
M]) . which satisfy the commutation relations 



[Qi, Q2] = Q2, [Q3, Q4] = 



(16) 



and thus the corresponding symmetry algebra is a real- 
ization of the algebra 2^2. According to the results of 
Patera and Winternitz 52| again, an optimal system of 
one-dimensional subalgebras is those generated by 

Q2, Qs, Q4, Qi + aQz, Qi + aQ4, Q2 + aQ4, Q2 + aQz. 

The corresponding similarity variables and reduced 
ODEs are listed in Table HH 

Table II. Reduced ODEs for /(«) = e^" 



(let A 7^ 0, E = (e^ 



"yyyyy)y)- 



i 


Subalgebra 


Ansatz u = 


y 


Reduced ODE^ 


5 


(Q2> 


■"(y) 


X 


E = 


6 


(Q3> 


v(y) + In X 


t 


144e^" -Xvy =0 


7 


(Q4> 


■"{y) 


t 


Vy=0 


8 


(Qi+aQ-i) 


v{y) + ^\nt 


xt-" 


E = - ayVy 


9 


{Qi+aQ4) 


v{y)-\\nt 


X -I- In t 


E = avy- 


10 


{Q2+aQ4) 


■"(y) 


X — at 


E = —aVy 


11 


{Q2+aQa) 


<y) - ^ 




E = ayVy - ^ 



C. f{u) = u"' (m / 0) 

In this case, the symmetry operators are Eq. (|15p . 
These operators share the same commutation relations 
Eq. ([T6l) . Hence an optimal system of one-dimensional 
subalgebras is the same as the case f{u) — e'^". The cor- 
responding similarity variables and reduced ODEs are 
listed in Table Hill 

Table III. Reduced ODEs for f{u) = u™ 

(let mjtO,E= (v"'Vyyyyy)y). 



i 


Subalgebra 


Ansatz u = 


y 


Reduced ODE^ 


12 


(Q2> 




■"(y) 


X 


E = 


13 

14 
15 
16 
17 
18 


(Q3> 

(Q4) 

(Qi +aQ3) 
{Qi+aQi) 
(Q2 +aQ4> 
(Q2 +aQ3> 


v{y)x'^ 

■"(y) 

6a-l 

i;(y)i ™ 
v{y)t~'^ 

■"{y) 

v{y)e~~^ 


t 
t 

xt-'^ 
X + alnt 
X — at 


+ ^ n {^+k)=Vy 

k = -A 

Vy=0 

E=^v-ayvy 
E = avy- 
E = —aVy 
E = ayVy - 



IV. INVARIANT SOLUTIONS 

Using the above reduced ODEs, we can construct some 
invariant solutions for the original equations (j4l). It is 
easy to see that some of the similarity variables in the 
tables U In] and IIIIl have a clear physical interpretation. 
Besides, for some higher order reduced ODEs, they can 
be further reduced by using new symmetries. Below, we 
discuss some facts related with some types of similarity 
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solutions with physical interest and obtain some partic- 
ular solutions. Different types of solutions are separately 
analyzed. 



A. Source and Sink Solutions 

There are two ODEs, i.e., ODEjjgjand ODI!|y^ among 
our reduced equations are related to this type of solu- 
tions. In fact, if we choose a = -p^^ in 0DE|y5], then the 
similarity solution has the form 



u{t, x) 



Thus, if m > —6 it is clear that u{t,x) — S{x) as t — ?> 
and the similarity solution is a source solution; if m < 
—6 it is clear that u{t,x) — > S{x) as i — > -l-oo and the 
similarity solution is a sink solution. Furthermore, we 
can also observe that, for the above choice of a = —tt;, 
ODEIj^ can be integrated once to obtain 



1 



V Vy 



m + 6 



yv = k, 



where k is an arbitrary constant. Thus, we can obtain a 
class of source solutions and sink solutions for the general 
thin- film equations (jj]) with f{u) = u™ (i.e., equation 
^) by solving the above fifth-order ODE. If we further 
choose k as zero, we have 



1 



^yyyyy 



m + 6 



yv 



This equation admits the symmetry group corresponding 
to the infinitesimal generator v — ydy + -^vdy. Taking 
into account that the invariants of its first prolongation 
and setting 



6 



Xl 



vy 



Ul 



^y..[yv' v) ^ 



this equation becomes a fourth-order ODE: 

- m^uluio-^xixixi + 5m'^{3muixi + 2mu\ 

- &ul)u\uix,x^x^ + IQm^ulul^^^^ 

- b'w?\p.lrn?u\^_^ -\- 20m{m — 3)M^uia;i 

+ {Irr? - 48m -(- 12)u\\uyUxxvxi + lOSm^uf^^ 
150m''(m - ■i)u\u\,^^ + \hnC'{lm^ - 48rn 

-f- 12)u\u\^^ \^m^{m - 3)(5to^ - 48m 
12)u\uY.^^ -f- [m5a;r'"+V(™ + 6) 

4- 72(m - 2)(m - 3)(m - 6)(2to - 3)xi]w? 

+ 12m{2m^ ~ SOm^ + 315m'^ 

- 720to -I- 540)m^ = 0. 

Thus, source and sink solutions can be also obtained by 
solving the above fourth-order ODE. 



If we choose m = —6 in ODEjYg], then the similarity 
solution has the form 

u(t,x) — rw(— — r), 

thus, if a > it is clear that u{t,x) — S{x) as i — )■ +oo 
and the similarity solution is a sink solution; if a < it 
is clear that u(t, x) — ^ S{x) as t — )■ —oo and the similarity 
solution is a source solution. As in the previous case, for 
the choice of m = —6, ODl!|j^ can be integrated once to 
obtain a fifth-order ODE 



v-^v. 



yyyyy 



ayv 



where k is an arbitrary constant. Consequently, source 
and sink solutions can be computed by solving a fifth- 
order ODE. 



B. Travelling- wave Solutions 

This types of solution corresponds to the reductions SI 
[TOlandfTTl In fact, in these three reductions the similarity 
variables are given by y = x — at, u — v, so that u(t, x) = 
v{x — at), thus the corresponding solutions are travUing- 
wave solutions. Due to the physical interest of this type 
of solutions, in what follows we study further symmetries 
of the associated ODEs and then construct some such 
kinds of solutions. First of all, we integrate these three 
equations once trivially and obtain 



OBE^:f{v)u 



vvvvv 



OD 



av = k, 
av — ki, 
av = k2, 



where k,ki,k2 are arbitrary constants. Because ODrTjQi 
and ODI'Iyyi are given by ODEjji for f{u) = e^' 
f{u) = u"* respectively, so we will focus on the OD^ 
below. This equation is invariant under the group 
translations in the y-direction, with infinitesimal genera- 
tor Set 
oy 



Xl =V, Ul =Vy , 

then the equation becomes a fourth-order ODE: 

\-^u\ux.j,^uxx^xvxx - uluix^xixtxi) + axiuf = kul- 
We further suppose that 

V = ^{xi,ui)dx^ + cj){xi,ui)dm 



(17) 



is an infinitesimal generator of the last equation, then the 
coefficients ^(a:i,ui) and (j){xi,ui) are satisfied with 



^ = axi + b, 

[5(a + c)axi + ba — (4a + 5c)k]f{xi) 
[—aaxi + [ak — ba)xi + kb]f'(xi) = 0. 



(18) 
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If fin) 



from the above system we can obtain 



5^ ~ (p — 0, which means that ODEmri has no nontrivial 
symmetry. Thus, it can not be reduced again. Conse- 
quently, the travelhng wave solutions for Eq. Q with 
f{u) — e'^" can be computed by solving a fourth-order 
ODE: 

l^u\ui^^ui^^xtxt - ufui^^^^^^^J + axiul = ku\. 
If f{u) = u™, then we have three cases from system 

m 



(i) k^O, C = 0, = 0; 

(m) fc = 0, m = 1, = axi + b, (p ~ — |aiti; 
(m) fc = 0, TO^l, i = axi, (f)-"^- 



5 -aui. 



Due to the triviality, the first case is excluded from the 
consideration. From the second case, we can get a ratio- 
nal travelling wave solution for Eq. Q with f{u) = u in 
the form 



i=0 



For the third case, we can set 



1^ 5 — m._, 

X2 — UiX-^ " , U2 — Xj^" [XiUi H — Ui) , 

o 

then Eq. (IT71) can be reduced to: 

625a;2M2"2a;2a:2a;2 " 125 [50x2^2:^2 -|- (11m 

- 25)a;2U2 + 75u2]a:2U2U2x22:2 + ^'^^bxlul^^ 
+ 125[3a;2U2(llm - 25) + 275]a:^U2uL2 

+ 25[125(5to - 12)a;2U2 + (46m2 - 225m 

+ 250)a;^M^ + 2%2h\x-2.ulu2x., + (2^ + 875to2 

- 250m3 - 1250m -I- 625a2:^ -I- 625)a;^u^ 
+ 10(48m3 - 375m2 + 875m - 625)a;^u^ 

125(38m2 - 195TO 22h)x\u\ + 13125(2m 

- 5)a;2U2 + 65625u^ = 0. 

Consequently, the travelling wave solutions for Eq. ^ 
with /(m) — u™'[m ^ 1) can be computed by solving a 
third-order ODE. 

Finally, we consider a special situation when f{u) ~ 
ue~" and k — 0, in which system (|18p infers that v — 
-hdxi + uidu^. Let 

X2=Mie~, M2 = -(5uixi -I- u)~ e ^, 
then Eq. ([TTll is reduced to: 

xlulu2x^X2X2 - xlu2{lOX2U2x2 + llulx2 

+ l5u2)u2x2X2 + 15a;2"L2 + lla;2"2(5 

-|- 3x2U2)ul^^ + X2u\{A6x\u\ + 125x2li2 



105)u2x2 + x^{'62^axl + 24)u^ + 96x^w^ 

t;2U2 



Therefore, the travelling wave solutions for equation ([4]) 
with f{u) = ue~" can be computed by solving a third- 
order ODE too. 



C. Waiting-time Solutions 

0DEIY3] is a first-order equation that can be easily 
solved, m this way we obtain a family of waiting-time 
solutions for the sixth-order thin film equation Q cor- 
responding to f{u) = M™ (if m ^ 3/2,2,3 or 6). These 
solutions are given by 



u(i, x) 



[m(<o-i) IT (— + ^)] ^^0, 

k=-4 

0, a; < 0. 



where Iq being an arbitrary constant. 



D. Blow-up Solutions 



ODEjg] is also a first-order equation. Solving it we get 
for the sixth-order thin film equation (U) with f{u) = e"^" 
the corresponding similarity solution 



^ 1, {X - Xpf 

^ ' ' A 144(to - t) 

where to is an arbitrary constant. This solution describes 
a localized blow-up at x — xq. Note that the solution is 



only valid if 



144(to-t) 



< 1, then, it ceases before t = tg. 



190x^1*^ + 210x2^2 + 105it2 = 



V. CONCLUDING REMARKS 

We have carried out a detailed group-theoretical anal- 
ysis for the generalized one-dimensional sixth-order thin 
film equation Q which arises in considering the motion 
of a thin film of viscous fiuid driven by an overlying elas- 
tic plate. A complete Lie point symmetry group clas- 
sification for the class ([4]) have been performed under 
the continuous equivalence transformation group. Based 
on these, a complete list of symmetry reductions of the 
classification cases have been derived by making use of 
the optimal system of one-dimensional subalgebras of the 
corresponding Lie symmetry algebras. Furthermore, in- 
variant solutions of the Eq. ^ with different functional 
form of / have been constructed by solving the reduced 
ODEs. In particular, by focusing our attention in those 
aspects with physical interest, we have found: 

1. The thin film equation (|4]), for the case f{u) — 
(which corresponds to equation ©), m > —6 ad- 
mits source solutions and m < —6 admits sink so- 
lutions. These solutions are related to the solutions 
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of a fourth-order ODE. If m = —6, Eq. (|4]) admits 
source and sink solutions. In this case these fami- 
hes of solutions are related to a fifth-order ODE. 

2. The thin film equation (|4]) has travelling- wave solu- 
tions. In the case f{u) = u™, for m = 1 the equa- 
tion admits a rational travelling- wave solutions, for 
m 7^ 1 the problem of finding these solutions can 
be transformed into the problem of solving third- 
order ODEs. In the case f{u) = e^", the travel- 
ling wave solutions can be computed by solving a 
fourth-order ODE. While for the case f{u) = ue^", 
the travelling wave solutions of equation (|4|) can be 
computed by solving a third-order ODE. 

3. Waiting-time solutions in the case f{u) = u™, 
and blow-up solutions in the case f{u) = e^" 
are obtained in the context of symmetry reduc- 
tions. However, it should be noted that these 
two types of solutions can also be obtained by 
means of variable separation. In the first case one 
takes u{t,x) — T{t)X{x) and in the second case 



u{t,x) = T{t)+X{x). 

These results may lead to further applications in physics 
and engineering such as tests in numerical solutions of 
Eq. Q and as trial functions for application of variational 
approach in the analysis of different perturbed versions 
of Eq. Q . Other topics including nonclassical symmetry, 
non-Lie exact solutions and physical applications of class 
(|4]) will be studied in subsequent publication. 
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